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Phase separation of colloid-polymer mixtures, described by the Asakura-Oosawa 
(AO) model, confined in a thin slit pore is studied by grand-canonical Monte Carlo 
simulation. While one wall is a hard wall for both particles, at the other wall there 
acts an additional repulsive potential on the colloids only. Varying the strength of 
this potential, a crossover from capillary condensation to an interface localization 
transition is found. The latter occurs under conditions where in the bulk almost 
complete phase separation has occurred. 
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I. INTRODUCTION AND MOTIVATION 



The current paradigm to create nanoscopic devices has enhanced the interest in the 
changes of phase behavior and structure when fluids are confined to linear dimensions that 
are no longer very large compared to the particle sizes. E.g., for a fluid (or fluid mixture) 
in a nanoscopic slit pore an interesting interplay of finite size effects, wetting or drying, 
and capillary condensation or eva por ation-like behavior can be expected 
R y. y. Iioi [ ll. 12, 13, 14, 15, 16, 17, A particularly interesting theoretical prediction 
[3j, |4j, y, ly, iTj concerns systems that can undergo phase separation (a fluid separating into 
vapor and liquid, or a binary mixture with a miscibility gap) confined between "competing 
walls" . By competing walls we mean that one wall favors one of the coexisting phases and 
the other wall favors the other one (the generic model Q, 0, S, O, B| is a thin ferromagnetic 
Ising film with oppositely oriented surface fields at both surfaces). When one brings the 
system in a state where phase separation occurs in the bulk, the thin film also exhibits 



pnase separation, witn an mtei 

B fl m n, 12, 15, 0, Q 



,3. 



)hase separation, with an interface parallel to the walls in the center of the film 

However, sufficiently far off from bulk criticality then a 



transition occurs where the interfacegets localized at the wall(s). This transition can be of 
first order or of second order and in the latter case it was argued to belong; 

to the universality class of the two-dimensional (2d) Ising model. However, while interfaces 
parallel to the walls of thin films could be observed experimentally in thin films of polymer 
mixtures [3], no experimental observation of this novel type of transition has as yet been 
reported. Indeed, the theoretical models have always invoked perfect symmetry between the 
two coexisting phases in the bulk (spin reversal symmetry of the Ising magnet [H U, 0, Ifl Q] , 
or strictly symmetric polymer blends 0,0,0,0,0], which are hard to realize in nature). 

In the present work we show that this interface localization should also be observable in 
colloid-polymer mixtures, which are known to be model systems for the experimental study 
of phase separation j^] and interfacial fluctuations 3, 22]. The large size of the colloid 
particles (of the order of 1 /xm) should allow to prepare slit pores which are of the order of 10 
- 100 particle diameters wide, with strongly repulsive walls which are essentially perfectly 
flat on this mesoscopic scale, modeling thus a "hard wall" boundary. However, coating a 
wall with a layer of long endgrafted flexible polymers of the same chemical type as used 
in the colloid-polymer mixture, one could create an additional repulsive interaction to the 
colloids to cancel (partially or fully) the depletion attraction jbll of the colloids to the 
hard walls. Varying the grafting density in the resulting "polymer brush" 24] and/or the 
molecular weight of the grafted chains the range and strength of this short range repulsion 
between colloids and the wall could be fine-tuned, within some limits 



In the following we shall present model ca 
mixtures, the Asakura-Oosawa (AO) model 



culations of a simple model for colloid-polymer 
confined between two inequivalent walls: a 



hard wall on one side, and a hard wall plus short range repulsion acting only on the colloids 
on the other side. Varying the strength of this repulsion we demonstrate a crossover from 
capillary condensation-like behavior to interface localization transitions. 



II. MODEL AND SIMULATION TECHNIQUE 



The colloids are modeled as hard spheres with diameter a c , the polymers are (soft) spheres 
with diameter a p . The polymers may not overlap with colloids, but there is no interaction 
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between the polymers. Recall that flexible polymers in solution form random coils with a 
rather large gyration radius which may interpenetrate at very low energy cost. We use a size 
ratio q = <j p /<j c = 0.8, since the phase diagram of this model in the bulk has been carefully 
studied previously [23, We choose a c = 1.0 as our unit of length. When there are 
N p polymers and N c colloids in the considered volume V, the polymer and colloid packing 
fractions are defined as r] p = na^N^/ {QV) and r/ c = ira^N c /(6V), respectively. The volume 
is given by V = LxLxD. Here, D is the width of the slit pore in ^-direction. Periodic 
boundary conditions are used in x- and ^/-directions parallel to the walls. The Lx L wall at 
z = simply is a hard wall, while the wall at z = D is described by an additional repulsive 
potential. It has the following form (we absorb a factor ksT here, /cs=Boltzmann's constant, 
T=temperature) : 

00 for z < 

e for 0<z< a c /2 (1) 
otherwise. 

The strength e of this wall-colloid repulsion is varied between and 2.5, while for polymers 
both walls are purely repulsive. Slit widths D = 5 and D = 10 were studied, while L was 
varied from L = 10 to L = 30, in the framework of a finite size scaling [2^ analysis. 

The simulation is carried out in the grand-canonical ensemble, with the chemical poten- 
tials /i p , /i c of polymers and colloids (or fugacities z p , z c ) as independent variables. Following 
common practice, we use the so-called "polymer reservoir packing fraction" r/ p = Tiz^a^/Q 
rather than z p as the temperature-like variable. As in the study of bulk critical behavior 
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2cS] and "capillary condensation" of the colloids on hard walls we use a. grand 



canonical cluster move |2a] together with a very efficient reweighting scheme, successive 
umbrella sampling [30], to obtain the distribution function PL(r] c \n T , z c ). This distribution 
function is defined as the probability to observe the system with colloid packing fraction r\ c 
at "inverse temperature" rf v and colloid fugacity r] c . For states far away from phase coexis- 
tence, Pl is a single-peaked function, while near phase coexistence a double-peak structure 
develops 2J Q- 2oL [30I . The precise location of the value of z c at which two-phase 

coexistence occurs is given by the equal weight rule 01 . The positions of the two peaks 
of Pl then yield (preliminary) estimates for the two branches of the coexistence curve, the 
"liquid" branch rf c and the "vapor" branch r\ v c . However, near criticality these estimates are 



affected by finite size effects (27, 



281 129[ . To deal with the latter, we study reduced moments 
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FIG. 1: Colloid concentration profiles r] c (z) and polymer concentration profiles r] p (z) as function of 
z for a thin film with asymmetric walls (hard wall at z = 0, while for the other wall at z = D = 10 
the potential of Eq. (JTJ acts, with e = 2.5). Profiles were obtained at tj c = 0.18, rf = 0.70 (a), 
r/ c = 0.18, rf p = 0.95 (b), rj c = 0.05, ^ = 1.20 (c), and r) c = 0.33, rf p = 1.20 (d). For profiles (c) and 
(d), the choices r] c = 0.05,0.33 roughly correspond to the two branches of the coexistence curve, 
see Fig. [51 

of Pl at phase coexistence, denning an analogue of the order parameter of the Ising model, 

oo 

m = r] c - (77 c ), (77c) = J r] c P L (r] c \r]l, z c )drf c , (2) 

o 

and higher order moments, (m p ) = J mP Pj^d^p, z c )dr] c . Following the behavior of ratios 
such as U4 = (m 2 ) 2 / (m 4 ) along the path in the (z c , 7/ p ) phase along which phase coexistence 
occurs for several choices of L one can estimate the critical point ("cumulant intersection 
method" 291). 



5 



1.3 
1.2 

_a 1.1 

1 .OH 

0.9 
0.8 




a) D=5, L=30 



0.0 



0.1 



0.2 



0.3 



0.4 




• e 


= 0.5 


n e 


= 1.5 


♦ e 


= 2.0 


4 e 


= 2.5 



0.4 



FIG. 2: Coexistence curves for D = 5, L = 30 (a) and D = 10, L = 20 (b), showing four choices of 
e in each case. The data points with open symbols show the vapor-like (7/") and liquid-like (77^) 
branches of the coexistence curve, as extracted from the peaks of Pl. The broken curves show the 
coexistence diameter 5 = (rj^ + n^)/2, on which the critical points are marked with full symbols. 
The critical values 7?p jCrit of rf v are highlighted by dotted horizontal lines. 

III. SIMULATION RESULTS 

Figure [T] shows typical concentration profiles for a slit of width D = 10 at r\ c = 0.18 and 
e = 2.5 and three choices of the polymer reservoir packing fraction: 7f p = 0.7 (in the one 
phase region of the bulk), rf p = 0.95 (in the two phase region of the bulk) and at rf p = 1.20 (as 
we shall see, this is in the two phase region of the thin film: therefore two profiles are shown 
here, corresponding to the two coexisting phases). For rf v = 0.7 one can see that there is an 
enhancement of the colloid concentration on the hard wall, as expected from the depletion 



attraction already noted in previous studies 



13. 



Mill 



On the repulsive wall, the colloid 



concentration is somewhat depressed (and the polymer concentration slightly enhanced), but 
in the center of the film both concentrations are roughly constant, as expected for bulk-like 
behavior. 

For rf p = 0.95, on the other hand, the profiles are very different: phase separation in a 

colloid-rich and a colloid pure phase has occurred, separated by an interface in the center of 

the slit. The interfacial profile resembles that of an interface between bulk coexisting phases 

fl 

(broadened by capillary waves [32]). For = 1.2, finally, the interface is localized either at 
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FIG. 3: Fourth-order cumulant U4 (a) and susceptibility (b) ksTx = L 2 D((m 2 ) — (|m|) 2 ) plotted 
vs. r]L, for a slit of thickness D = 5 and e = 1.0, for several values of L as indicated in the figure. 
Insert in a) shows the slope Y\ of the cumulants, evaluated at 77p Crit = 0.937 ± 0.003, on log-log 



scales to extract the exponent v (note 
the extrapolation of the peak positions, from which if crit 
(a). 



2£| Y\ oc L l l u ), yielding \jv ~ 1.056. The inset in b) shows 
0.935 ± 0.005 results, consistent with 




1r 



FIG. 4: Phase boundary » mf T (D, e) plotted vs. rf p for several choices of e. For comparison, also 
the symmetric slit, e = pjj] and the bulk phase boundary fi^^ 1^2 ] are included. Full symbols 

Mcoex 



mark the locations of the critical points. Inset shows A/i = //|?oex ~~ Mcoex (D,e) versus rf . 
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FIG. 5: Dependence of the critical value r/ p 

crit on e i f° r D — 5 (circles) and D — 10 (squares). 

the right wall or at the left wall (Fig. [U;,d): phase separation has occurred in the thin film 
(Fig.H. 

Since Fig. |21 refers to finite L, the data for the coexistence curves do not merge at the 
critical points, but rather extend somewhat into the one-phase region ("finite size tails"). 
However, repeating the simulation for several choices of L and performing a finite size scaling 
I29I ] a naly sis, as done for this model both in the bulk (2?! Issl and in thin films with symmetric 
walls [l7l ll8 | . the critical point can be determined rather reliably. Figure El gives an example 
for e = 1.0. Due to crossover effects the cumulant intersections do not occur precisely in a 
point, but rather are spread out over some region, but nevertheless, the critical point r/p Crit 
can be inferred with reasonable accuracy. The effective exponents v c r pa 0.947 (extracted 
from Y 1 , Fig. I2ta)) and (/3/V)cff ~ 0.155 (extracted from A = (rjf. — rj^)/2 at r/p crit {not 
shown}) deviate somewhat from the 2d Ising values {y — 1, /3/v — 0.125), similarly as in 
thin Ising films with competing walls 

When e increases, rf p crit is shifted to considerably larger values and at the same time the 
chemical potential /x coex (L>,e) moves towards the bulk curve, /i^oex (Fig- EI) , and is almost 



reached for e = 2.5. The limiting value ^p jCrit (-D — > 00) 
would yield an estimate for the wetting transition ry pw 
roughly, we estimate rf p vf pa 1.20 ± 0.05. 



'or which u rn ^(D — ► 00, e) = [A 



bulk 

coex 



12j of this model. Very 



Comparing to the case of symmetric hard walls (e = 0), a non-monotonic variation of 
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?7p Crit on e is found (Figs. EJ)- A similar behavior has been seen for asymmetric polymer 
mixtures studied by self-consistent field theory ^(|. Our finding implies that this non- 
monotonic variation is a rather general phenomenon, and it holds beyond mean field theory. 

IV. CONCLUSIONS 

In a colloid-polymer mixture confined by hard walls, depletion attraction leads to the 
formation of colloid-rich layers near the walls. By adding a suitable short range interaction 
at one of the walls, polymers can be effectively attracted to this wall whereas the depletion 
attraction of colloids is reduced at it. In this case of a colloid-polymer mixture confined 
between asymmetric walls, our simulations demonstrate that an interface localization transi- 
tion can occur. Our results also imply that semi-infinite colloid-polymer mixtures, confined 
by a hard wall, should exhibit a wetting transition at sufficiently large polymer fugacity. We 
suggest that interface localization transitions could be realized experimentally by suitable 
coating of one wall of a slit pore with a polymer brush. 
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